Operational determination of multi-qubit entanglement classes via tuning of local 

operations 
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We present a physical setup with which it is possible to produce arbitrary symmetric long-lived 
multiqubit entangled states in the internal ground levels of photon emitters, including the paradig- 
matic GHZ and W states. In the case of three emitters, where each tripartite entangled state belongs 
to one of two well-defined entanglement classes, we prove a one-to-one correspondence between well- 
defined sets of experimental parameters, i.e., locally tunable polarizer orientations, and multiqubit 
entanglement classes inside the symmetric subspace. 
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Entanglement is a distinctive property of quantum 
physics associated with the nonseparable character of 
multipartite quantum systems. For the case of two-qubit 
systems, entanglement is well understood and can be pre- 
cisely quantified Apart from the trivial disentangled 
case, three qubits possess two genuine tripartite inequiv- 
alent entanglement classes 0, 0] • Efforts have been done 
recently towards higher number of qubits 0, [H, 0], in- 
cluding an inductive method though so far no com- 
prehensive and scalable classification has been developed. 
In this letter, we introduce a physical setup consisting of 
N emitters, incoherently radiating single photons that 
may be absorbed remotely by detectors equipped with 
polarizers and producing long-lived multiqubit entangled 
states among the emitters. We show that it is possible 
to associate well-defined sets of locally tuned polarizer 
orientations with multiqubit entanglement classes, allow- 
ing their monitoring in an operational manner. Hereby, 
multipath quantum interferences, associated with qubit 
permutation symmetry, play a key role in explaining the 
underlying physics. 

We consider a chain of N equally separated single pho- 
ton emitters, say trapped neutral atoms, trapped ions, 
quantum dots, or any other equivalent physical system 
with access to similar behaviour. Each emitter defines a 
three-level A system, where \e) denotes the excited state 
and the two long-lived sublevels, |+) and |— ), define a 
qubit. We assume that the transitions between the ex- 
cited state and the two lower sublevels have an equal 
wavenumber and dipole moment, and that they are cir- 
cularly polarized, cr + and cr_, respectively. Figure 1 ex- 
emplifies the Af-emitter case discussed throughout this 
paper. All emitters are initially excited and we will study 
the cases in which all spontaneously emitted photons 
are detected by N detectors located in the far-field re- 
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FIG. 1: (color online). Proposed experimental arrangement. 
N excited emitters are aligned in a row, each of them defin- 
ing a three-level A-system. A long-lived entangled state is 
obtained in the N emitter qubits after detecting the N spon- 
taneously emitted photons with N detectors equipped with 
polarizers. The final A-qubit state is tuned and determined 
by the polarizer orientations. 



gion, each of them being equipped with a polarization 
filter in front. The far-field detection ensures the era- 
sure of which- way information of the arriving photons, 
and the polarizers allow the generation of quantum su- 
perpositions of the lower atomic states when considering 
arbitrary polarizations. As a consequence each photode- 
tection event projects the emitters onto linear combina- 
tions of the long-lived states |+) and | — ) [8]. This results 
at the end in a coherent superposition between the qubit 
states |zb, . . . , ±). The indeterminacy of which detector 
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state of the W type ^ , the weakly entangled state 



FIG. 2: (color online). Pyramid of entanglement paths 
for the case of 3 emitters initialized in the excited state 
|e,e, e). The figure illustrates the intermediate states, hor- 
izontally displayed, during three successive photodetection 
events realized by three detectors di, cfe, and cfe, equipped 
with general elliptical polarizers oriented along ei, €2 and 
€3 (ei = (Xi<J+ + f3i<T-), respectively. After each detection, 
the unnormalized global state of the system is the sum of all 
states in the circles weighed by the tagged prefactors. Left- 
down arrows denote <x + transitions, and right-down arrows 
represent <r_ ones. The final state components are shown 
inside colored frames and the colored arrows represent the 
quantum paths leading to these different components. Only 
the red circle states are obtained via a single quantum path, 
while the blue and green ones are the result of three different 
interfering quantum paths. 



has projected which emitter in its ground states implies 
the existence of many quantum pathways between the ini- 
tial fully excited atomic state and each of the final state 
components. This produces a multipath quantum inter- 
ference effect [9] that we will tune by modifying the polar- 
izer orientations. We remark that several experimental 
setups may be in condition to implement the concepts 
introduced in this paper [Io|, EH, EH EH . 

As will be shown explicitly below, it is always possi- 
ble to find suitable polarizer orientations to produce any 
desired state totally symmetric with respect to permuta- 
tions of the emitters. Hereby, linear polarizers allow the 
generation of a state of the GHZ type [lj] , the maximally 
entangled state 
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(|l ,O ,...,O ) + ... + |O ,...,O ,l )), 

(3) 



with |O ) = (|+)-e^|-))/V / 2. In the 3-qubit case, where 
each tripartite entangled state belongs to one of two well- 
defined entanglement classes [2|, [3j : GHZ or W family, 
it turns out that even the rotation of a single polarizer 
allows us to switch from one entanglement class prepa- 
ration to the other, which are univocally determined by 
the number of distinct polarizer orientations in the ex- 
perimental setting. Those polarizer manipulations can 
be said to be local with respect to the polarizer /detector 
positions, though they are not local with respect to the 
qubit positions: it is well known that one cannot con- 
vert states from different families into each others using 
stochastic local operations and classical communication 
(SLOCC) [15[. Remarkably, the above described local 
polarization rotations will make possible the transitions 



S class <-> W class <-> GHZ class, 



(4) 



where "S" stems from "separable". 

There are several physical systems where the genera- 
tion of GHZ and W states with three or more qubits have 
been experimentally achieved: in trapped ions \wi Ell, 
in Rydberg atoms crossing microwave cavities [18|, and 
in photonic systems [EH, [2(j. Furthermore, other mul- 
tiqubit entangled states have been realized in different 

with different pur- 
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physical setups 
poses [26] and potential applications in quantum informa- 
tion. Though several paradigmatic entangled states have 
been produced in the lab, there is no study, to our knowl- 
edge, that associates operationally given experimental 
configurations with multipartite entanglement classes. 

The outlined behavior can be understood from the ex- 
plicit calculation of the different states of the TV-emitter 
system after the successive photon detection events. A 
properly located detector in the far-field with a gen- 
eral elliptical polarizer oriented in the xy plane of the 
circularly polarized light along the polarization vector 
e = acr + + (3<j- with arbitrary complex coefficients a 



and (3 ( 
tion [27 



1) can implement the operator ac- 



\GRZ N ) = -= (|+, ...,+) + e»\- ...,-», (1) 
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V (e) = a J2\+) j (e\+ l 3j2\-)j(e\, 
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with arbitrary relative phase (/>, or a separable (product) 
state 



ISiv) 



(2) 



with 1 10 ) = (|+)+e^| — ))/v // 2. An intermediate polarizer 
configuration permits the generation of the multiqubit 



up to an insignificant prefactor. Here, the sum over j 
runs over all emitters and |±)j(e| is the projection oper- 
ator from state \e) to state |±) for emitter j. Starting 
with N emitters in their excited state |e, . . . , e), the de- 
tection of the N emitted photons by N detectors with 
polarizer configuration €1, . . . , (e^ = c^cr+ + A^ 7 -) 
projects the emitter system onto the final state \ipf) = 
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FIG. 3: (color online). Operational monitoring of tripartite 
entanglement classes by changing the polarization configura- 
tions. If the three polarizer orientations are all different, the 
final state belongs to the GHZ class; if two polarizer orienta- 
tions are different, the final state belongs to the W class; if 
all polarizers are identically oriented, the final state belongs 
to the S class. 



)/V2 oriented along angles 
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2iV + (fc " 1) iV' 
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or any configuration resulting from permuting the in- 
dices. Here, the term inside the square brackets is only 
present for the case of an even number of emitters. The 
state |GHZtv) appears naturally when using all distinct 
polarizer orientations uniformly shared out over 2tt. 

In contrast, we are left with the separable state |Sjv) 
of Eq. ([2]) with linear polarizers all identically oriented 
along the angle 0/2 : 



V(€n) • • • P(ei)|e, . . . , e). Due to the symmetry proper- 
ties of those operators, the final state as well as all inter- 
mediate states are totally symmetric with respect to per- 
mutations of the emitters. If entanglement is produced at 
the end of the detection process, we are ensured that only 
genuine multipartite entangled states belonging to the 
accessible symmetric entanglement classes will be gen- 
erated. The intermediate states produced in the course 
of the successive photon detection events can be ordered 
in a pyramidal manner, displaying the various possible 
quantum paths towards the generation of the desired fi- 
nal state (see Fig. [2] which exemplifies the case N = 3). 
At the end of the N photon detection process, the sys- 
tem is found in a coherent superposition of all possible 
product states of the 7V-qubit system in the |±) basis. 
All probability amplitudes related to the different quan- 
tum paths add up and yield interference terms (tagged 
prefactors on Fig. [2]). The final state \ipf) is found to be 
a linear combination of all symmetric Dicke states with 
k |— ) excitations [HiSH, \D N (k)) : 



(9) 



N 



(6) 



where AT is a normalization prefactor, and 

{C k N f /2 E Pi 1 -.-Pi k a ik+1 ...a iN , (7) 



Ck 



with C^- the binomial coefficient of N and k. Thereby, 
we can generate any symmetric state with respect to per- 
mutations of the emitters using suitable elliptical polar- 
izer orientations. This is a direct application of Vieta's 
formulas [30[ : an arbitrary symmetrical state can be ex- 
panded in the Dicke state basis as Y^k dk\D]\[(k)) and is 
produced in our setup using K polarizers oriented along 
vectors with ai/fa identifying to the K roots of the 

polynomial P(z) = £f (-l) K ~ k ^C%/C§d k z k with K 
the degree of this polynomial, the remaining polarizers 
oriented along <r + . 

For example, the maximally entangled state | GHZ n) 
of Eq. (pQ) is generated with linear polarizers 7rg k = 



The state | Wat) is also generated with linear polarizers, 
all except one oriented identically along <j>/2 with the last 
orthogonal to the N — 1 first, 



= . . . = On- 



Oat = — - zb — . 
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(10) 



or any configuration obtained with permutation of the 
indices. 

For three emitters, those particular results, associating 
the states | S3) , IW3), and IGHZ3) with certain choices of 
polarization angles, suggest a wider physical picture. In 
fact, we can associate in an operational manner specific 
polarizer configurations of our proposed physical setup 
with the three paradigmatic entanglement classes ap- 
pearing in the 3-qubit case : the number of distinct po- 
larizer orientations identifies univocally the entanglement 
class, be S, W, or GHZ, as shown in Fig. [3j 

According to Diir et al. [2] , the GHZ class is formed by 
states characterized by a non vanishing 3-tangle [3l| , the 
W class by states with a zero 3-tangle and non-zero single 
qubit von Neumann entropies, while the separable state 
class is characterized by zero values of these entanglement 
measures. In our proposed scheme, the 3-tangle r of the 
final state \ipf) reads 

r = ^AfVi/i - a 2 pi\ 2 \a 1 p 3 ~ 0^3/^1 1 2 |o^ 2 /3 3 - a 3 f3 2 \ 2 

(11) 

and vanishes only when 2 polarizers are equally oriented. 
In this case, the local entropies vanish only when the 
third polarizer coincides with the two first. When switch- 
ing from a configuration with three distinct polarizer ori- 
entations to a configuration with three identical ones, via 
the intermediate case where two of them are equal, we 
transit successively from the GHZ class, to the W class, 
and end in the S class, as seen in Fig. [3j This remarkably 
shows the potentiality of the proposed setup for associat- 
ing operationally a physical setting with 3-qubit entan- 
glement SLOCC classes of states symmetric with respect 
to permutations of the emitters. These results encourage 
a possible generalization to the arbitrary Af-qubit state. 
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Finally, we give an estimation of the expected fidelity of 
our scheme, e.g., for generating the |GHZ 4 ) state using 4 
adjacent ions localized in a linear trap. Assuming a 5 jam 
separation between ions with a 5 nm confinement in the 
transverse direction and an azimuthal detection window 
of 1°, we estimate a fidelity of about 90% for the genera- 
tion of the 4-qubit state, largely above the 50% minimal 
value for proving genuine four-qubit entanglement of the 
state [32]. With usual ellipsometry techniques, polarizer 
orientation uncertainties can be made insignificant com- 
pared to the finite size effect of the detection window. 
With an excitation rate of several tens of MHz [8], the 
counting rate of the needed fourfold coincident events 
can be of the order of several tenths of Hz using CCD 
cameras covering a fair area in the detection plane. In 
general, the counting rate decreases with the number of 
qubits. This might limit the scalability of the scheme as 
in similar experiments related to entanglement produc- 
tion @, ei, Ei, im . 

Note added .-After submission of this work, we became 
aware of a nice experimental observation of an entire fam- 
ily of four-photon entangled states [331 ]. 

L.L. thanks Alexander von Humboldt Foundation for 
the funding through a Humboldt Research Fellowship. 
G.S.A. thanks NSF grant NSF-CCF-0524673 for sup- 
porting this collaboration. C.T. and J.v.Z. gratefully ac- 
knowledge financial support by the Staedtler foundation. 

E. S. thanks to Ikerbasque Foundation, EU EuroSQIP 
project, and UPV-EHU Grant GIU07/40. T.B. thanks 

F. Bastin, P. Mathonet, P. Lecomte and M. Rigo from 
University of Liege for helpful discussions. 



[1] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 
[2] W. Diir, G. Vidal, and J. I. Cirac, Phys. Rev. A 62, 
062314 (2000). 

[3] A. Acm, D. Brufi, M. Lewenstein, and A. Sanpera, Phys. 

Rev. Lett. 87, 040401 (2001). 
[4] F. Verstraete, J. Dehaene, B. De Moor, and H. Ver- 

schelde, Phys. Rev. A 65, 052112 (2002). 
[5] L. Lamata, J. Leon, D. Salgado, and E. Solano, Phys. 

Rev. A 75, 022318 (2007). 
[6] L. Chen and Y.-X. Chen, Phys. Rev. A 74, 062310 (2006). 
[7] L. Lamata, J. Leon, D. Salgado, and E. Solano, Phys. 

Rev. A 74, 052336 (2006). 
[8] B. B. Blinov, D. L. Moehring, L.-M. Duan, and C. Mon- 
roe, Nature 428, 153 (2004). 
[9] C. Cabrillo et a/., Phys. Rev. A 59, 1025 (1999); X.-L. 

Feng et al, Phys. Rev. Lett. 90, 217902 (2003); C.-S. Yu 

et a/., Phys. Rev. A 75, 044301 (2007). 
[10] U. Eichmann et al, Phys. Rev. Lett. 70, 2359 (1993). 
[11] Y. Miroshnychenko et a/., Nature 442, 151 (2006). 
[12] J. Beugnon, M. P. A. Jones, J. Dingjan, B. Darquie, G. 

Messin, A. Browaeys, and P. Grangier, Nature 440, 779 



(2006). 

[13] D. L. Moehring, P. Maunz, S. Olmschenk, K. C. Younge, 
D. N. Matsukevich, L.-M. Duan, and C. Monroe, Nature 
449, 68 (2007). 

D. M. Greenberger, M. Home, A. Shimony, and A. 
Zeilinger, Am. J. Phys. 58, 1131 (1990). 
Only approximated conversions can be realized. See e. 
g. P. Walther, K. J. Resch, and A. Zeilinger, Phys. Rev. 
Lett. 94, 240501 (2005). 
D. Leibfried et al, Nature 438, 639 (2005). 
H. Haffner et a/., Nature 438, 643 (2005). 

A. Rauschenbeutel, G. Nogues, S. Osnaghi, P. Bertet, M. 
Brune, J.-M. Raymond, and S. Haroche, Science 288, 
2024 (2000). 

J.-W. Pan, D. Bouwmeester, M. Daniell, H. Weinfurter, 
and A. Zeilinger, Nature 403, 515 (2000). 

B. P. Lanyon and N. K. Langford, quant-ph/0802.3161 
(2008). 

J.-W. Pan, M. Daniell, S. Gasparoni, G. Weihs, and A. 
Zeilinger, Phys. Rev. Lett. 86, 4435 (2001). 
Z. Zhao, Y.-A. Chen, A.-N. Zhang, T. Yang, H. J. 
Briegel, and J.-W. Pan, Nature 430, 54 (2004). 

C. A. Sackett et a/., Nature 404, 256 (2000). 

D. Leibfried, M. D. Barrett, T. Schatz, J. Britton, J. 
Chiaverini, W. M. Itano, J. D. Jost, C. Langer, and D. 
J. Wineland, Science 304, 1476 (2004). 
N. Kiesel, C. Schmid, G. Toth, E. Solano, and H. Wein- 
furter, Phys. Rev. Lett. 98, 063604 (2007). 
P. Kok, H. Lee, and J. P. Dowling, Phys. Rev. A 65, 
052104 (2002). 

G. S. Agarwal, J. von Zanthier, C. Skornia, and H. 
Walther, Phys. Rev. A 65, 053826 (2002). 
In an N spin-^ system, the Dicke states, usually denoted 
by \S,m), are defined as the simultaneous eigenstates of 

- 2 

both the square of the total spin operator S and its 
^-component S z , where S(S + l)h 2 and mh are the cor- 
responding eigenvalues. The N + 1 states with the high- 
est value S = N/2 form a special subset of Dicke states 
that are symmetric under qubit permutation, and can 
be characterized by the number of excitations k and de- 
noted \DN(k)). For details see R. Dicke, Phys. Rev. 93, 
99 (1954). 

[29] C. Thiel, J. von Zanthier, T. Bastin, E. Solano, and G. 

S. Agarwal, Phys. Rev. Lett. 99, 193602 (2007). 
[30] Vieta's formulas are sum-of-roots and product- 

of-roots identities for polynomial. See e. g. 

E. W. Weisstein, "Vieta's Formulas." From 
MathWorld - A Wolfram Web Resource. 



[14 
[15 



[16 
[17 
[18 



[19 

[20; 

[21 

[22; 

[23' 
[24 

[25; 

[26; 

[27 
[28 



\http://mathworld. wolfram.com/V ietasFormulas. html 



[31] The 3-tangle is a measure of genuine tripartite entangle- 
ment that is invariant under permutations. For a formal 
definition see V. Coffman, J. Kundu, and W. K. Woot- 
ters, Phys. Rev. A 61, 052306 (2000). 
[32] M. Seevinck and J. Uffink, Phys. Rev. A 65, 012107 
(2001). 

[33] W. Wieczorek, C. Schmid, N. Kiesel, R. Pohlner, O. 
Giihne, and H. Weinfurter, Phys. Rev. Lett. 101, 010503 
(2008). 



